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1. INTRODUCTION 

The design of controller for mobile robots are studied and developed in many works [1]—[3]. Fuzzy 
controllers are designed using evolutionary algorithms and fuzzy approach [4]-[6]. Fuzzy controllers have 
been successfully used for practical systems [7]-[10]. It is suitable for controlling systems without need for 
mathematical models. Several evolutionary algorithms such as genetic algorithm (GA), artificialbee colony 
(ABC), and particle swarm optimization (PSO) are adopted for the design of an optimal fuzzy controllers. 
GA was widely used for fuzzy controller design and applied in different applications [11], [12]. PSO is also 
used for fuzzy controller design and applied in different applications [13], [14]. ABC is utilized for fuzzy 
controller design and applied in different applications [15], [16]. It was applied for numerical function 
approximation in [17]. Zhu and Kwong [18] proposed a chaotic bee colony algorithm. A modified ABC was 
presented to enhance the exploitation phase [19], [20]. A comparison of ABC with PSO was presented in 
[21]. A modified ABC methods based on PSO were presented in [22]—[24]. 

In this article a new improved ABC is presented to be used in the design of an efficient fuzzy 
controller for motion tracking of mobile robot. It determines appropriate membership functions. The 
proposed improved ABC is based on works presented in [23]—[25]. It uses some features of the PSO 
algorithm to improve the performance. The PSO algorithm determine the best result through adjusting the 
speed and location of each individual. The modified ABC (MABC) make a balance between the exploration 
and the exploitation which results in going through the global space and increase the speed of convergence. 
The modified algorithm also use a chaotic initial population and an extra random search is used to improve 
the best results for the current iteration in the scout phase. The rest of this paper is arranged as follows: in 
section 2, mathematical model is presented. In section 3 the artificial bee colony and the proposed modified 
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ABC are presented. In section 4 a modified ABC based fuzzy control scheme design method is explained. 
Results are presented in section 5 and section 6 is for conclusion. 


2. THE SUGGESTED METHOD 
2.1. Mobile robot 
In general mobile robots are described by two models: 
a. Kinematics model 
Kinematic model is used for small speed, small load, and small acceleration. The linear velocity v 
and the angular velocity Ê are expressed by the following [5]: 


ll = f A H a) 


in (1), r is the wheel radius, distance between wheels is 2b, 0 represents the direction angle. wr is the right 
wheel speed, wl is the leftwheel speed. 
b. Dynamic model 

Figure 1 shows the mobile robot diagram. It consists of two wheels. Both are placed on one axis. 


The front wheel is not actuated and it is used for balancing. The state vector is defined as x = [v 0 òf, the 
output vector is y = [v 0 ]” and the manipulated variable is u = [u, u,]’. The state-space equations are [5]: 


x = Ax + Bu (2) 
y = Cx (3) 
A, B, and C are: 

a, 0 0 b, by 
A=|0 0 |=| o |.c= i i p| with, 

0 0 ag b> —b, 

-2c 2 22d? o L kr b= krl 

ay = (Mswr2421,,)’ 42 ~ (yr2+2Iwl2)? 1 (Ms«r24+2Ilw)’ 2 (lyr2+2Iwl2) 


where 1, is the robot moment of inertia, M is the mass of robot, c is the distance between wheels, I, is the 
momemt of ineria of wheels, c is the viscous friction and r is tha radius of wheel. 


Figure 1. Mobile robot model 


2.2. Modification of artificial bee colony algorithm 
2.2.1. The ABC algorithm 

In ABC bees are in three groups: employed bees, onlookers, and scouts. The work of the employed 
bee is the search for available food sources. The original artificial bee colony concept can be described by the 
foraging behaviors of bee colony. If there exist two food places: A and B. At the starting, a bee forager will 
begin as an employed bee. The bee with no knowledge about places of food has two options [17]: 
— It works as a scout and begins searching for a place of food. 
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— It works as a recruit when seeing the dances and search for a food place. If it finds a food, it store the 
location in memory and begins exploiting it, and it becomes as an employed forager. It takes a quantity of 
nectar and returns back to put the nectar in a store. Then the bee has to select one of the following: 

— It may become as a nonaligned follower. 

— It may dance to recruit other bees. 

— It may keep to forage at the food place. 

Mathematically, the original ABC algorithm may be stated as follows. In the beginning the ABC 
algorithm creates random initial food locations of N solutions, where N is number of employedbees. The 
solution x; (i= 1, 2, ..., N) is M x1 vector. Then the nectar amount fitness(fiti) for each individual is found. In 
this algorithm, the value of function fitness is the nectar quantity. In the employedbees phase, the new found 
food source i near the current food x; is determined by (4) [16], [17]: 


Zij = Xi + Oi Xiyy — Xx) (4) 


where both k and j are random integers k E (1, 2, N) and j E (1, 2, M) and k is not equal i, Ø;j is a random 
value in the range [-1,1]. The new generated individual is compared with the current one to keep the better 
using a greedy choice method. In the onlookerbees phase, the food source is selected a ccording to a 
probability which depends on the fitness of afood shared by employedbees. Probability is determined by: 


P= =f (5) 


N fit; 
Xj=1fitj 


In the scoutbee phase, if there is no improvement in a food source when reach the limit, it is deleted and the 
employed bee may work as a scout. The scoutbee generates a new extra random food location using: 


H =x) in +rand[0,1] (xax — XÍ in) (9 


x min min 


min and x Jax are smaller and larger values of x /, respectively. The steps are repeated for the 
assumed number of repetitions, or until a stop condition is reached. 


where x/ 


2.2.2. Modified ABC 

The modified ABC algorithm used here is based on the methods proposed in [24], [25] for 
overcoming the disadvantages of the original algorithm. The first modification is the use of some features of 
the PSO algorithm to enhance the solution quality. In the employedbee phase, speeds and locations of 
individuals are updated by (7) and (9), respectively similar to that in PSO. 


vi (k +1) = wevi(k) +c «17, *[ PPOs! — x) (k)] + cz * r * [4° — x) (k)] (7) 


The weight w modify the particle momentum by changing the weight for the former speed [25], the 
value of w is proposed here to be updated via an exponentially increasing function as: 


w= (Wmax = Wpmin)e ~ (Mmariter-iter)/( iter-1)) a 
and 
xi (k +1) = xik) + v (k +1) is 


In the case of minimization, the particle of personal best solution at the step after is given by: 


pre (k)iff@(k+D)= f (pP (k+1)) 


best k+1)= 10 
Pr ED) 5 det) if F Oy (R+D) < f (pI (k+ 0) aes 

gbest can be found from (11): 
ge = min {f (p ?®)} ,i=1,2,...,N, (11) 
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The second modification in the MABC algorithm is to create an initial population by chaotic dynamics. The 


chaotic approach is used to enhance the diversity of the population and increase the convergence speed. The 
logistic map is used to create a chaotic value as: 


ch(i + 1) =4ch() * (1-ch(i)) (12) 
Then convert ch to x, a chaotic value in the duration [x; x,] by (13): 

xG) =x, + chy — x) (13) 
xı and x, are the smaller and larger value of x, respectively. The third modification is to add an extra random 
search in scout phase of the modified algorithm. S-random numbers in the interval [xl xu] are generated to 
enhance the search activity and to avoid falling in local optimum: 

Rd(i) = x, + rand[0,1](x, — x,),i=1,..S (14) 
where S is the length of the random sequence. Finally, use the next equation to calculate a new solution x: 

XG) = (1 — a)gbest + a Rd(@i),i = 1,...... S (15) 
where æ is the decaying factor defined here as: 


a= e(- (iter—1)/(maxiter—iter)) (16) 


It can be noticed that a decreases with the increase of iterations. 


The Modified ABC algorithm 


1. Select the population size N individual vector length M, maxiter, Wmax, Wmin, Vmax, aNd Vmin. 
2. Find a chaotic initial population using (13) {xi; i=l, . . N}. Determine f {fi; i = 1,. 
N}. Find gbest and pbest i. 

3. While stop is not satisfied do 

semployed phase 

4. Fori= 1 to M do 


5. Update velocities and positions using (7) and (9). Find the new best personal 
solution p?@*(k+1) using (10) and the global best solution using (11) 

6. Endfor 

sonlookerphase 


7. Fori = 1 to N do 

8. If rand < Prob (i) 

9. Update speed and position using (7) and (9). Find the best personal solution using (10 
and update trail. 

0. Endif 

1. If traili exceed limit, then 

% scoutphase 

2. change xi to the new chaotic value from (13). 

3. Endif 

4. Find the gbest. 

5. Begin searching using Random process S times using (14) and from (15) find the new 
candidate soultion %. Redetermine the gbest. 

6. iter = iter +1 

7. Endwhile 


3. MODIFIED ABC BASED FUZZY CONTROL OF MOBILE ROBOT 
In this method an optimization steps of the fuzzy controller was made by MABC. Each particle 
represents the complete fuzzy system parameters. The following points are assumed: 

— All values are normalized to lie between -4 and 4. 

— Apexes of first and last membership functions (MFS) are at —4 and 4 respectively. 

— In the middle triangular membership functions are used and, in the edges, trapezoidal membership 
functions are used. 

— The number of MFS is an odd integer (seven MFS). They are symmetrical about vertical axis. With the 
above assumptions the design parameters are the centers of MFS. The number of parameters to be 
optimized for each input or output variable is three (symmetrical MFS about the y-axis). The aim is to 
find torques u; and u, such that the two-wheeled mobile robot can follow the required path. The proposed 
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fuzzy control scheme shown in Figure 2 consists of two 2-input- l-output controllers named velocity 
controller with input the velocity error e, and the change in error Ae, and angle controller with inputs the 
azimuth error eg and Aeg where: 

ey = Va — V and €g = 8g — 9 
Ae (t) =e,(t) —e, (t — 1) and Aeg(t) = eg(t) — eg (t— 1) 


The output of these two fuzzy controller y1 and y2 are used to find torques u, and uj: 


Ur = yı + Y2 ; U = y1 — Y2 


e Speed Controller ; 
Angle Controller : 


Mobile Robot 


Figure 2. Block diagram for motion control 


where vq , Oq are the required speed and the required angle, respectively, v, 0 are the actual speed and the 
actual angle of the robot, respectively. Seven values: negative big (NB), negative middle (NM), negative 
small (NS), zero (ZO), positive small (PS), positive middle (PM), and positive big (PB) are used for each 
variable. Outputs of both controllers are computed using the weighted average method [3]. The MABC 
algorithm is applied for optimization of the two fuzzy controllers. The population consists of N particles and 
expressed by: 


P= [p1p2 .......pN | 


Membership functions are symmetrical so that only three centers values are required thus, three 
groups of three center values [all, a12, al3] and [a21, a22, a23] and [b11, b12, b13] for the two input 
variables and the output variable of the velocity controller and [a31, a32, a33] and [a41, a42, a43] and 
[b21, b22, b23] for the variable of the angle controller. The vector of particle in the MABC algorithm 
consists 18 elements, 9 elements for each fuzzy controller. The vector of the particle at any generation can be 
expressed by: 


P = [p1p2 p3 p4 p5 p6 p7 p8 p9 p10 p11 p12 p13 p14 p15 p16 p17 p18} 
= [a11 a12 a13 a21 a22 a23 b11 b12 b13 a31 a32 a33 a41 a42 a43 b21 b22 b23] 


and the particle h at generation g is denoted by Pp ; 


4. EXPERIMENTAL RESULTS 

The aim is to test the performance of the developed controlscheme. The mobile robot shown in 
Figure 1 is considered. It has three wheels, two active wheels and one wheel for balance. Parameters values 
are from [5] and shown in Table 1. 


Table 1. Values of mobile robot parameters 
Parameters Values __ Units 


l, 10 kg.m? 
M 200 kg 
l 0.3 m 
ly 0.005 kg.m? 
c 0.05 kg/s 
r 0.1 m 
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To simulate the above mobile robot model, MATLAB and Simulink is adopted. The search space 
consists of 18 dimensions, nine dimensions specified for velocity controller and another nine dimensions for 
azimuth controller. The MABC initial parameters values are given in Table 2 [24]. The mobile robot is tested 
for tracking circular path to study the way in which the proposed control scheme steers the mobile robot to 
track a circular path with desired velocity described by vg = 0.3 [meter/sec] and desired azimuth given by: 
Oa = (2 * 3.14 * f(t)/m) rad, with m=-5 where m is the slop, f(t)=t ,0<t<5 and T,= 0.01 sec. In the test, the 
mobile robot is assumed to start at [Xo yo 90]=[0.3 m 0 m 0 rad]. Results are presented. Table 3 shows the rule 
base. Figure 3 shows the optimized membership functions for the two input variables and the output of the 
controller using MABC. The membership functions are symmetrical about y-axis. Figure 4 shows the mean 
square error against iterations for original ABC and MABC. Figure 5 shows the circular trajectory. Figure 6 
and Figure 7 show error in motion x(t) and y(t) respectively. Clearly results show good improvement when 
MABC is used. 


Table 2. MABC initial values 


Value Parameters 
50 Size: no of bees 

80 Maximum iterationnumber 
18 Dimension 

1.4 Parameter c; 

1.4 Parameter c2 

0.9 Wmax 

0.4 Wmin 

-1, 1 Vmin, Vmax 


Table 3. Rule base for the fuzzy controller 
eAe NB NM _ NS ZO PS PM PB 
NB NB NB NB NB NM NS ZO 


Membership Function Plots 
à NB NM NS z PS PM PB 
0.5 
-4 -3 -2 -1 input vardbie “error 1 2 3 4 
Membership function plots 
, NB NM NS z PS PM PB 
0.5 
= a a! oO 
4 3 2 inpul variable “Change a Error” 2 3 4 
Membership function plots 
NGS NM NB z PS PM PB 
t a a m a . 
j 
] } 
os | 
| | i; | 
a 4 a E a d 
4 -3 -2 2 3 4 


-4 C] 4 
output variable "output 


Figure 3. MFS using MABC 
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10? 
—— Using ABC 
| —— Using MABC 


fitness values 


iterations 


Figure 4. Mean square error versus iterations 


—— desired trajectory 


Figure 5. The circular trajectory 
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Figure 6. Error in motion X(t) 
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Figure 7. Error in motion Y(t) 


CONCLUSION 
This paper has presented an optimal motion fuzzy control scheme using the proposed modified ABC 


to achieve trajectory tracking of mobile robot. The modified ABC gives good results. The tracking 
performance is better when the modified ABC is adopted. The additional good facilities added to the original 
ABC algorithm make the modified algorithm more powerful. The tracking mean square error is much 
reduced with the modified algorithm. 
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